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ABSTRACT 

We propose a new scenario for supermassive star (SMS; > 10^ Mq) formation in 
shocked regions of colliding cold accretion flows near the centers of first galaxies. 
Recent numerical simulations indicate that assembly of a typical first galaxy with 
virial temperature Tvir > 10'' K proceeds via cold and dense flows penetrating deep 
to the center, where the supersonic streams collide each other to develop a hot (~ 10* 
K) and dense (~ 10"^ cm~'^) shocked gas. The post-shock layer first cools by efhcient 
Lya emission and contracts isobarically until ~ 8000 K. Whether the layer continues 
the isobaric contraction depends on the density at this moment: if the density is 
high enough for coUisionally exciting H2 rovibrational levels (> 10* cm~^), enhanced 
H2 coUisional dissociation suppresses the gas to cool further. In this case, the layer 
fragments into massive (> 10^ Mq) clouds, which collapse isothermally (~ 8000 K) 
by the Lya cooling without subsequent fragmentation. As an outcome, SMSs are 
expected to form and evolve eventually to seeds of supermassive black holes (SMBH). 
By calculating thermal evolution of the post-shock gas, we delimit the range of post- 
shock conditions for the SMS formation, which can be expressed as: T > 6000 K 
(nn/lO* cm-3)-i f^^. < ^^-3 g^^^ > qqqq _ qqqq ^ > ^^-3^ 

depending somewhat on initial ionization degree. We found that metal enrichment 
does not affect the above condition for metallicity below ~ 10"^ Zq if metals are in 
the gas phase, while condensation of several percent of metals into dust decreases this 
critical value of metallicity by an order of magnitude. Unlike the previously proposed 
scenario for SMS formation, which postulates extremely strong ultraviolet radiation 
to quench H2 cooling, our scenario here naturally explains the SMBH seed formation 
in the assembly process of the first galaxies, even without such a strong radiation. 

Key words: stars: formation. Population III — dark ages, reionization, first stars 
— galaxies: formation, nuclei 



1 INTRODUCTION 

Discovery of high-2 quasars has demonstrated the existence 
of supermassive black holes (SMBH) with ~ 10^ Mq al- 
ready at the age of the universe < 1 Gyr (e.g.. Fan 2006; 
Willott et al. 2007; Mortlock et al. 2011). As an origin of 
such SMBHs, seed BH formation as a remnant of population 
III stars (Msecd ~ 100 Mq) and their subsequent growth by 
merger and gas accretion have been studied by a number 
of authors (e.g., Haiman & Loeb 2001; Volonteri, Haardt 
& Madau 2003; Li et al. 2007). Even with the Eddington 
accretion rate MEdd = ^Edd/ec^, where Z/Edd is the Edding- 
ton luminosity and e ~ 0.1 is the radiative efficiency, the 
growth time tgrow ~ 0.05 ln(A/BH/Msccd) Gyr becomes as 
long as the age of the universe at 2: ~ 6 (~ 0.8 Gyr): the 
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seed BHs are thus required to keep growing at the Eddington 
rate. However, negative feedbacks by the growing BHs pre- 
vent such efficient accretion (Milosavljevic, Couch & Bromm 
2009; Alvarez, Wise & Abel 2009; Jeon et al. 2011). 

As a solution to this, alternative possibility of massive 
seed BH formation by direct collapse of supermassive stars 
(SMS; > 10^ ^0) li^-s been considered by some authors. 
Specifically, SMS formation in massive halos (Tvir > 10* K) 
irradiated with strong far ultraviolet (FUV) radiation has 
often been studied (e.g., Bromm & Loeb 2003; Regan & 
Haehnelt 2009a, b; Shang, Bryan, & Haiman 2010). Since 
H2 molecule, the main coolant in the primordial gas, is pho- 
todissociated with strong FUV radiation in the Lyman and 
Werner bands, clouds under such an environment collapse 
isothermally at ~ 8000 K by Lya cooling without fragmenta- 
tion, if they are massive enough with > 10* Mq . As an out- 
come of such collapse, SMSs are expected to form. A massive 
seed BH as a remnant of SMS collapse reduces the growth 
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time to 10^ Mq within 0.46 Gyr and mitigates the growth- 
time problem by a big margin. This scenario, however, has a 
serious drawback: for this mechanism of SMS formation to 
work, extremely strong FUV radiation J^j^ > 10^ — 10'^ (in a 
unit of 10~^^ erg s^^ cm^^ Hz~^ sr~^) is required (Omukai 
2001; Bromm & Loeb 2003; Shang, Bryan & Haiman 2010), 
while the fraction of halos irradiated with such intense FUV 
fields with J^^ > 10^ is estimated to be < 10"'' at 2 ~ 10 
(Dijkstra et al. 2008), i.e., only extremely rare halos satisfy 
the condition for SMS formation. Moreover, if high energy 
components, such as cosmic rays or X-rays, are present along 
with the FUV radiation, their ionization effect promotes the 
H2 formation and then strongly suppresses the SMS forma- 
tion (Inayoshi & Omukai 2011). Although the above scenario 
might be still viable considering the rarity of high-z SMBHs, 
it is worthwhile to explore another possibility. 



In this paper, we propose a new scenario for SMS for- 
mation in post-shock gas of cold accretion flows in form- 
ing first galaxies. Recent numerical simulations of galaxy 
formation have revealed that, in halos with virial tempera- 
ture Tvir > 10* K, the shock position does not stay at the 
virial radius and shrinks inside owing to the efficient Lya 
cooling, and the accreting cold gas penetrates deep to the 
center through dense filamentary flows (Birnboim & Dekel 
2003; Keres et al. 2005; Dekel & Birnboim 2006; Dekel et al. 
2009; Bromm & Yoshida 2011). The supersonic flows col- 
lide each other and the resultant shock develops a hot and 
dense (~ 10* K and ~ 10"^ cm~"^) gas near the center (Wise 
& Abel 2007; Greif et al. 2008; Wise, Turk & Abel 2008). 
By studying thermal evolution of the shocked gas, we have 
found that, if the post-shock density is high enough for the 
H2 rovibrational levels to reach the local thermodynamic 
equilibrium (LTE), the efficient coUisional dissociation sup- 
presses H2 cooling, and the gas cannot cool below several 
thousand K. Massive clouds with > 10''' Mq formed by 
fragmentation of the post-shock layer subsequently collapse 
isothermally at ~ 8000 K by the Lya cooling. Without fur- 
ther fragmentation, monolithic collapse of the clouds results 
in the SMS formation. Note that, unlike the previous SMS 
formation mechanism, strong FUV radiation is not required 
in our scenario. Similar analysis has also been carried out 
by Safranek-Shrader, Bromm & Milosavljevic (2010), who 
studied the fragmentation of the cold-stream shocked layer 
considering the effects of radiation field and chemical enrich- 
ment, but for a single post-shock condition (4 x lO^cm""^, 
1.1 X 10* K, and the equilibrium chemical abundances). 



The organization of this paper is as follows. In Section 
2, we describe the model for calculation of thermal evolution 
in the shocked gas. In Section 3, we present our results and 
clarify the conditions for isothermal collapse leading to the 
SMS formation in terms of post-shock density and tempera- 
ture. Effects of metal enrichment is also considered here. In 
Section 4, we analyze thermal processes in the shocked gas 
and discuss the reason for the bifurcation of thermal evo- 
lution in more detail. Finally, we summarize our study and 
present some discussions in Section 5. 



2 MODEL 

In this section, we describe our model for calculation of ther- 
mal evolution in hot and dense shocked regions formed by 
collision of cold accretion flows in the first galaxies. 

2.1 Evolution in the post-shock layer 

We consider the thermal evolution in the post-shock layer 
under the assumption that the flow is steady and plane- 
parallel. Since the post-shock temperature is as high as the 
virial temperature of first-galaxy forming halos (Tvir > 10* 
K), cooling by the Lya emission is efficient early on. The 
post-shock flow is compressed almost isobarically as long as 
the gas cools effectively (e.g., Shapiro & Kang 1987). Within 
the steady-state approximation, the conservation of mass 
and momentum leads to the following relationships between 
the density po, pressure po, and flow velocity vq just behind 
the shock with those in the post-shock flow p, p, and v: 

pv = povo, (1) 

pv^ + p = povo + po. (2) 

Along this flow, we solve the energy equation 
dE p + E dp 

where E is the internal energy per unit volume, d/dt is the 
Lagrangian time derivative, and Anct is the net cooling rate 
per unit volume. Assuming a strong shock and neglecting 
the thermal pressure in the pre-shock flow, po — Spo^'o 
is satisfled just behind the shock front. Thus, we approx- 
imate the right-hand side of equation ((2]) by 4po^'o- The 
cooling term Anct includes the radiative cooling by H, II2 
and HD, and cooling/heating associated with chemical re- 
actions. We solve the chemical reactions of primordial gas 
among the following 14 species; H, II2, e~, H"'", Hj, H~, 
He, He+, He++, D, HD, D+, HD+, and D". We adopt the 
same coefficients for the cooling/heating and the chemical 
reactions as in Inayoshi & Omukai (2011) except for omit- 
ting the radiative/cosmic-ray ionization and dissociation in 
this calculation. In studying effects of metal enrichment, we 
add the cooling by the fine-structure-line emission of Cn 
and Oi to the primordial processes described above. As- 
suming the fraction of metals depleted to dust grains to 
be the same as in the Galactic interstellar gas, we set the 
number fractions of C and O nuclei in the gas phase with 
respect to H nuclei to a;c,gas = 0.927 x 10~*(.Z/Zq) and 
2:0, gas = 3.568 X 10~*(Z/Zq) (Pollack et al. 1994). We fol- 
low HoUenbach & McKee (1989) in calculating the cooling 
rates of Cn and Oi. We curtail the C and O chemistry by 
simply assuming that all the C and O are in the states of Cn 
and Oi, respectively, from the following consideration: with 
lower ionization energy (11.26 eV) than H atom, C is pho- 
toionized by weak background radiation and in the state of 
Cii, while O is in ionization equilibrium with H and almost 
neutral for < 8000 K, where the Oi cooling is important. 
Molecular cooling of metals (e.g., CO and H2O) is not in- 
cluded since its cooling is not important in the temperature 
range relevant for the bifurcation of thermal evolution (> 
several 10^ K). 

Next, we consider the condition for the gravitational in- 
stability during the isobaric compression of the post-shock 
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layer and thus for its fragmentation. For the isobaric com- 
pression, the dynamical time tdyn(= p/ (dp/dt)), which char- 
acterizes thermal evolution, is approximately equal to the 
cooling time 



(3/2)nHfcBr 

Anct 



(4) 



where wh is the number density of H nuclei and T is the 
temperature. On the other hand, the growth timescale for 
the gravitational instability is given by the free-fall time 
(e.g., Larson 1985) 



32 
37rGp' 



(5) 



As long as the cooling is effective enough and so tcooi(— 
^dyn) ^ ta, the post-shock layer continues to be compressed 
isobarically. Ifowever, once the cooling becomes ineffective 
and tcooi exceeds ta, the contraction of the layer halts and a 
dense layer begins to develop inside the post-shock region. 
For example, since the growth rate of the baryonic mass 
in halos of ~ 10* Mq at z ~ 10 is ~ 4 x 10"^ MQ/yr 
(Dekel et al. 2009), a gas of > 10'' Mq can accumulate in 
> 3 X 10^ yrs. If sufficient gas supply is available through 
the accretion flow, the dense layer eventually satisfies the 
Jeans criterion, i.e., the sound-crossing time fcross becomes 
longer than the free-fall time (fcross iff)- The layer then 
fragments by the gravitational instability to produce Jeans- 
scale clouds. In this paper, assuming that the continuous 
gas supply to induce the Jeans instability is available, we 
estimate the fragmentation mass scale by the condition ta ~ 
icooi in accordance with Yamada & Nishi (1998). 



2.2 Evolution after fragmentation 

After fragmentation, the cloud collapses with its self-gravity, 
and its evolution cannot be modelled as a steady fiow any- 
more. Density evolution in a cloud collapsing by the self- 
gravity is described by the Penston-Larson self-similar solu- 
tion (Penston 1969; Larson 1969), which has density profile 
with a flat core of the Jeans scale and an envelope with the 
power- law density distribution p(r) oc r~^. The density in 
the central core roughly increases in the free-fall timescale. 
We here calculate the evolution in the central-core part by 
using a one-zone model, where the density evolution is given 
by 



dp 
'dt 



P_ 

ts' 



(6) 



Namely, after the condition for fragmentation (tcooi ^ t«) 
is satisfied, we switch the density evolution described by 
equations ([T])-((2]) to that by equation ((HJ in our calculation 
and solve equation ([3]) for this density evolution. 

When the collapse proceeds significantly and the cloud 
becomes optically thick, the radiative cooling becomes inef- 
fective due to photon trapping. We assume the radius of the 
core to be half a Jeans length 



Aj 

2 V Gppmu ' 



(7) 



where fj, is the mean molecular weight. Since we consider 
the core of the collapsing cloud, the column density of i-th 



species is given by Ni = xinuRc, where xi is its concentra- 
tion. Using this value, we estimate the optical depth and the 
reduction rate of radiative cooling as in Inayoshi & Omukai 
(2011). 



2.3 Initial conditions 

According to numerical simulations of the first galaxy for- 
mation (e.g., Greif et al. 2008; Wise et al. 2008), the pre- 
shock number density and temperature of cold fiows and 
the shock velocity are typically 10^ cm~^, 200 K, and 20 km 
respectively, which correspond to the post-shock density 



4 X 10^ cm~^ and temperature 9000 K. With those fiducial 
values in mind, we carry out calculations for a wide range of 
initial number density and temperature: 10^ cm"'^ < nu,o < 
10'' cm-^ and 3000 K < To < 10^ K. Since H2, the main 
coolant below 8000 K, forms through the electron-catalyzed 
reactions 



R + e ->H +7, 
H--fH^H2-he-, 



(8) 



the initial ionization degree Xb,o, along with the initial H2 
concentration xh2,o, is important quantities for the subse- 
quent thermal evolution. In reference to the results of Kang 
& Shapiro (1992), who studied the chemical abundances in 
the pre-shock gas considering photo-ionization and dissoci- 
ation by UV radiation emitted from the shock, we regard 
3:^0,0 ^ 10~^ and xh2,o ~ 10~^ as typical ionization degree 
and molecular fraction, respectively. However, since cold ac- 
cretion fiows are far denser (~ 10^ cm"'') than the range 
Kang & Shapiro (1992) assumed (< 10"^ cm"^), the elec- 
tron recombination as well as the shielding of the UV photo- 
ionization/dissociation probably lower the pre-shock ioniza- 
tion degree 2:0,0 and elevate molecular fraction XH2 from 
those values. Taking this uncertainty into account, we study 
the cases with a wide range of initial ionization degree and 
H2 fraction: 10'^ ^ a;c 10"^ and lO"'' ^ xh. ^ . 



3 RESULTS 

In this section, we present our results for thermal evolution 
after the gas experiences the cold accretion shock. We first 
consider the cases of primordial gas and then discuss effects 
of small metal enrichment. 



3.1 Primordial-gas case 

In Fig. [1] we show the temperature evolution of primor- 
dial gas for four post-shock conditions, indicated by two 
each open and filled circles. We here set the initial ion- 
ization degree and molecular fraction to a;o,o ~ 10~^ and 
x:h2,o = 10~®, respectively. First, we see the cases from the 
open-circle initial conditions in Fig. [1] whose temperature 
evolution is shown by the dashed lines. In the lower initial- 
density case (nH,o, ro) = (5 x 10^ cm~^, 3.6 x 10* K) among 
them, the post-shock gas cools by the Lya emission and is 
compressed isobarically. Although the Lya cooling becomes 
inefficient below 8000 K, enough H2 for cooling has already 
been formed by this time, which enables further tempera- 
ture decrease. HD is formed abundantly for < 150 K, whose 
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number density nH (cm-3) 



Figure 1. The temperature evolution of primordial gas after 
heated by the cold accretion shock for the initial ionization de- 
gree Xe,o = 10"'^ and H2 fraction x^^ q = 10~^. The evolution- 
ary tracks are shown with dashed and dash-dotted lines for four 
combinations of initial temperature and density, indicated by two 
each open and filled circles. From low-density or low-temperature 
initial conditions (dashed lines starting from the open circles), 
the temperature decreases below ~ 100 K owing to the H2 and 
HD cooling. On the other hand, from dense and hot initial condi- 
tions (dash-dotted lines from the filled circles), the clouds do not 
cool below ~ 8000 K and subsequently collapse almost isother- 
mally by the H atomic cooling. The triangle symbol on each track 
indicates the epoch when the post-shock layer fragments by the 
gravitational instability. The thick solid line, the domain above 
which is hatched, divides the initial conditions leading to these 
two different ways of thermal evolution. The two thin solid lines 
show the temperature evolution by the H atomic cooling (upper) 
and the H2 and HD cooling (lower) , respectively. These two evolu- 
tionary tracks are calculated by the density evolution of equation 
^ from the initial conditions Hh.o = 10 cm" -3, To = lO"* K, and 
Xefi = 2 X 10~* (upper) and lO"'^ (lower), respectively. For the 
upper track, the H2 cooling rate is set to zero. The diagonal dot- 
ted lines (lower-left to upper-right) indicate the constant Jeans 
masses, whose values are denoted by numbers in the Figure. 

cooling eventually lowers the temperature to ~ 50 K. At 
this point, without efficient coolant anymore, tcooi becomes 
longer than tg. Clouds with Jeans mass of several 10 Mq 
are produced by the gravitational instability. Also, in the 
case of lower initial temperature (10^ cm~"^, 3.3 x 10^ K), 
abundant H2 is formed immediately. The cooling by H2 and 
then by HD allows the temperature to plummet isobarically 
until ~ 100 K, where the fragmentation mass scale of a few 
10 Mq is imprinted. In both the cases, temperature evolu- 
tion after fragmentation converges to the well-known track 
for clouds collapsing by the self-gravity and cooling by H2 
and HD (lower thin solid line; Uehara & Inutsuka 2000; Na- 
gakura & Omukai 2005). 

Next, we see the two cases starting from the fiUed-circle 
initial conditions in Fig. [1] whose evolutionary tracks are 
indicated by the dash-dotted lines. In the higher initial tem- 
perature case (10* cm~"^, 3.6 x 10* K), the gas cools isobar- 
ically until 8000 K as in the open-circle cases. At ~ 8000 K, 
however, the density exceeds ~ 10* cm~'^, the critical den- 
sity for H2 to reach LTE. For higher density, H2 is rapidly 
dissociated coUisionally from the excited rovibrational lev- 
els, and thus sufficient H2 for cooling is never formed. Con- 



Figure 2. The effect of initial ionization degree on the range of 
post-shock initial conditions leading to isothermal collapse by H 
atomic cooling and thus SMS formation. The lines present their 
boundaries for the cases with different initial ionization degrees 
of Xe,o = 10-1 (soUd), 10-2 (long-dashed), IQ-^ (short-dashed), 
10^* (dotted), and 10~^ (dash-dotted), respectively. The molec- 
ular fraction ajjj., = 10-^ for all cases. 

sequently, the gas cannot cool below ~ 8000 K, and massive 
clouds with > lO"" Mq are formed by fragmentation. Also, in 
the case of initial temperature somewhat lower than 8000 K 
(2x 10^ cm-^, 5x10^ K), H2 cooling is suppressed by the col- 
lisional dissociation and fragmentation occurs immediately 
producing massive clouds with > 10^ whose tempera- 

ture increases to 8000 K by the compressional heating in the 
course of gravitational collapse. In both cases, the massive 
clouds thereafter collapse almost isothermally by the Lya 
cooling until very high density (~ 10^^ cm-^), where the 
cloud becomes optically thick to the H- bound-free absorp- 
tion (Omukai 2001). Such isothermally contracting clouds do 
not fragment in the later phase and thus collapse monolith- 
ically to SMSs, which eventually evolve to seeds of SMBHs 
(Bromm & Loeb 2003; Shang, Bryan & Haiman 2010). 

As seen above, the behaviors of thermal evolution can 
be classified into two types. The thick solid line in Fig. [T] 
corresponds to the boundary of initial conditions, from the 
above or below which subsequent thermal evolution bifur- 
cates. Namely, the post-shock conditions above the bound- 
ary (the hatched region) lead to the isothermal evolution 
at ~ 8000 K, while those below it result in the isobaric 
temperature decrease until < 100 K. The boundary on the 
low-density side {nu,o ^ 10* cm^'') can be fitted as 

rQ>8xio^( '^^f'" — 3)"'k. (9) 

\ 7 x 10^ cm--^ / 
This means that, after isobaric cooling to 8000 K, if the den- 
sity exceeds the H2 critical value for LTE (~ 10* cm-'^), the 
gas cannot continue further isobaric compression and starts 
isothermal collapse. For higher densities nH,o > 10* cm-^, 
the boundary is given by 

ro>5xio^( !^"-°_3 )"°'k. (10) 

V lO'' cm^/ 

In Section 4, we discuss physical processes determining the 
location of the boundary in more detail. 

We mention the effect of initial chemical composition 
on thermal evolution. In Fig. (2] we show the boundaries 
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for the SMS-forming conditions for different initial ioniza- 
tion degrees (10^''' ^ Xcfi ^ 10^^). The positions of the 
boundaries are almost independent of Xc,o for > 10* cm^'^, 
while the portions at < 10* cm~^ move to lower density 
with decreasing Xefi- In particular, for a^e.o = 10~^, this re- 
sults in the spiky domain around 8000 K extending as low as 
~ 3x 10^ cm~^. With the higher initial ionization degree, the 
more H2 is formed by the electron-catalyzed reactions ([8]), 
which results in the wider range of post-shock conditions for 
H2 cooling, i.e., the smaller range for SMS formation. The 
boundary for Xc,o ^ 10"^ asymptotically approaches that 
for higher a::c,o at > 3 x 10* K since the ionization degree 
jumps up immediately to ~ 10~^, even with smaller initial 
value, by effective coUisional ionization 

H-he" ^ H+ -f 2 e", (11) 

in this temperature range. We also studied the cases with 
different molecular fraction xh2,o ~ 10-^10-* and 10'^, 
and found that the boundaries for SMS formation is almost 
independent of xii2,o- This is because, even with high initial 
value, H2 is rapidly dissociated coUisionally for < 10*cm~'^, 
and its fraction reaches the equilibrium one, which is inde- 
pendent of initial value. 




number density nH (cm-3) 



Figure 3. The effect of metallicity on the range of post-shock 
initial conditions leading to SMS formation. The lines present 
their boundaries for the cases of metallicity Z = (solid), 
10-* (long-dashed), 5 X lO"" (short-dashed), 10"^ (dotted), 
2 X 10~^Zq (dot-dashed), respectively. With metallicity higher 
than Zcr = 10~^ Zq, the boundary remarkably shifts toward 
higher density. The initial ionization degree and molecular frac- 
tion are x^fi = 10~^ and xh^ q = 10~®, respectively. 



3.2 Metallicity effect 

Next, we consider the cases with slight metal enrichment. 
In Fig. [S] we present the boundaries for SMS-forming ini- 
tial conditions for various metallicities with ^ Z ^ 
2x10"^ Zq. With some metals, cooling by the fine-structure 
line emission of Cn and Oi can exceed that by H2 and plays 
an important role in thermal evolution. With metallicity as 
low as Z < 5 X 10~* Zq, the metal cooling does not affect 
thermal evolution around 8000 K and the boundary does 
not move from the primordial case. As seen in Section 3.1, 
from initial conditions hotter and denser than the boundary 
(i.e., the hatched region of Fig. [T]), enough H2 for cooling 
is not formed and only massive clouds are produced, which 
collapse isothermally thereafter. With increasing metallicity, 
the metal-line cooling becomes able to make the gas to cool 
below ~ 5000 K isobarically even without H2. Once the tem- 
perature decreases and the coUisional dissociation becomes 
ineffective, abundant H2 eventually forms and the gas cools 
to < 100 K by its cooling. By this additional cooling, the 
boundary of the SMS-forming initial conditions moves to 
higher density. With metallicity as high as Z ^ 10"'^ Zq, 
the cooling rate by Cn and Oi becomes comparable to the 
compressional heating at ~ 8000 K and ~ 10* cm~"^. There- 
fore, even without help of H2 cooling, the metal cooling alone 
is able to lower the temperature to the range where the H2 
coUisional dissociation is ineffective, and thus the boundary 
shifts to higher density. 

In summary, with metallicity higher than the critical 
value Zcr ~ 10~^ Zq, the boundary density becomes far 
higher than the typical post-shock value by the cold accre- 
tion shock ~ 10^ cm~^ and such an initial condition would 
be very difficult to be realized. Thus the possibility of SMS 
formation is strongly reduced for higher metallicity. On the 
other hand, as long as Z < Zcr, the range of initial condi- 
tions for SMS formation remains the same as in the primor- 
dial case. 



4 MECHANISM FOR THE BIFURCATION OF 
THERMAL EVOLUTION 

In this Section, we explain what processes are responsible 
for the bifurcation of thermal evolution in the region where 
the cold accretion shock is thermalized, and give a physical 
interpretation for the location of the bifurcation boundary 
for the post-shock conditions in Figs. [T] and [2] 

Efficient Lya cooling drives temperature in a hot gas 
rapidly to ~ 8000 K, where its cooling rate is sharply cut off 
as the atomic hydrogen is not excited for lower temperature. 
In Fig. |4l we show the cut-off temperature (dotted line), 
below which H2 takes over the role of the dominant coolant. 
Thus, for the post-shock gas to continue the isobaric cooling 
below 8000 K, the H2 cooling must become effective and keep 
the cooling time tcooi shorter than the free-fall time tg. The 
cooling time by the H2 cooling is given by 



tc 



(3/2)fcBr 



(12) 



where ~ AH2/("-Ha:^H2) is the cooling rate per an H2 



molecule (erg 



and has the density dependence 



(1 + nH.cr/nn)"^ (nn,, 
sity for LTE). In conditions under consideration (below the 
dotted line of Fig. |4} , the H2 fraction a;H2 is set by the equi- 
librium between the electron-catalyzed formation reaction 
and coUisional dissociation reaction Q and so 



^form 



(13) 



^ Only with high Xc ~ 10 *, the charge exchange reaction 
(H2 H+ ^ H2+ + H) becomes the main dissociation reaction of 
H2 at < 10^ cm""^. However, since the bifurcation boundary for 
^e,0 = 10~^ locates at higher density (> 10* cm~'^), the charge 
exchange reaction does not have any influences on the location of 
the boundary. Thus, we adopt equation I I13I I even for Xc ^ 10~*. 
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where fcform (H + e~ — >■ + 7) and feed (H2 + H — >• 3H) 
are reaction rate coefficients for the indicated reactions, re- 
spectively. Note that kcd depends on the fraction of H2 in 
excited states and thus on the density. The rate coefficient 
kcd significantly increases with density near the H2 critical 
density nn.cr, which results in the rapid decrease of by 
the effective coUisional dissociation at > 10^ cm~"^. 

Furthermore, since XH2 is proportional to x^, the re- 
combination also can be a key reaction to determine tcooi 
(cx x-'^) at < 8000 K. Since the recombination time tree 
(= l/arecnHXc; Oroc is the recombination rate coefficient) 
has the same dependence on Xe as tcooi, the ratio of these 
two timescales becomes independent of Xc and is approxi- 
mately given by 

icool _ (3/2)nHfcBr fccdQrcc 
tree 



0.9 



V5000 K 



104 cm- 



(14) 



The last expression above is valid for the density and tem- 
perature around tcooi — tree and > 10* cm~^, and the large 
temperature dependence of equation (|14p is due to that of 
feed. In Fig|4l we show the range of the parameters satisfying 
tcooi > tree (the hatched region), where the recombination 
effectively works during the isobaric contraction. Note that 
the hatched region appears only at > 10'^ cm~^, where k^d 
is enhanced significantly. On the other hand, the ionization 
degree is frozen during the isobaric compression at the den- 
sity lower than the hatched region as fcooi < tree- Under 
constant x^, the cooling time becomes shorter and shorter 
with decreasing temperature as t^ooi oc T®'^ (at ~ 5000 K) 
because the coUisional dissociation is strongly suppressed for 
lower temperature. Thus the efficient H2 cooling and resul- 
tant isobaric evolution continue until ~ 100 K. 

We first consider the cases with x^o 10~* and de- 
fer the discussion for the lower ionization cases later. As 
an example, in Fig. U we present the evolutionary tracks of 
the shocked gas with a;c,o = 10~^ for different initial tem- 
peratures (solid lines) and the bifurcation boundary of the 
SMS-forming initial condition (right dashed line). While the 
gas passes through the hatched region where t^cc > tcooi, 
the ionization degree and thus H2 fraction fall rapidly. If 
such a gas runs out of H2 before reaching the region where 
tcooi < tree, the condltiou for fragmentation teooi tg is 
immediately satisfied and the clouds formed in this way col- 
lapse isothermally thereafter (the solid lines starting from 
the filled circles in Fig|4]). However, there is a small mar- 
gin of the initial parameter range above the line tcooi ~ tree 
from which the gas manages to maintain a small fraction of 
H2 and can reach the region where teooi < tree- In this case, 
the post-shock layer can continue the isobaric contraction 
until ^ 100 K (the solid line starting from the open circle in 
FigHI- In summary, if the fragmentation condition tcooi ^ tft 
is met in the range where tree < tcooi, the post-shock layer 
cannot cool further and clouds formed at this moment be- 
gin the isothermal collapse. The set of initial conditions from 
which the evolutionary tracks meet the condition tg = icooi 
just on the line tree = tcooi corresponds to the boundary for 
the SMS formation in Fig.[T]in the high density regime (i.e., 
equation llUp . Note that the boundary given by equation 
l|10P reflects the density-temperature relation of tree ~ tcooi, 
which is mainly determined by the temperature dependence 




tcool=tff 

(xe,o=10'') 



10 10" 
number density nH (cm-3) 

Figure 4. The diagram showing relevant processes in setting 
the bifurcation of the thermal evolution of the shocked gas. The 
solid lines show the evolutionary tracks of the shocked gas with 



Xe,0 



10" 



whose post-shock conditions are nH,o 



2 X 10* 



cm-3 and To = 1.2 X 10*, 6.6 X 10^, and 5.1 X 10^ K, respectively. 
Two dashed lines present the bifurcation boundaries of the SMS- 
forming initial conditions for Xc,o = 10~^ and 10~^, respectively. 
The dotted line is the same as the upper thin line in Fig.[T] which 
corresponds to the cut-off temperature below which the H2 line- 
emission becomes the main cooling process instead of the hya 
emission. In the hatched region, the ionization degree rapidly falls 
by effective recombination within the cooling time of H2 (tcool > 
tree)- The dot-dashed line presents the condition of = t^ool 
for Xe = 10~^. Here, the cooling time tcool is evaluated by using 
equation 1131 1. 



of kcd- Due to the strong dependence of kcd on tempera- 
ture, the H2 is coUisionally dissociated very efficiently for 
the temperature higher than given by equation (|10|l . On the 
other hand, if the density is lower than the H2 critical den- 
sity nH,cr (~ 10* cm~^) after cooling isobarically to 8000 
K, the gas continues further isobaric contraction to several 
100 K by the ineffective coUisional dissociation of H2. This 
explains the fact that the low-density side of the boundary 
of SMS-forming parameters in Fig. [T] (i.e., equation [SJ cor- 
responds to the isobaric contraction track whose density at 
~ 8000 K is ~ 10* cm-^ 

Next, we consider the low ionization cases with Xc,o ^ 
10~*. As seen in Fig.|4l the portion of the bifurcation bound- 
ary (left solid line) for Xcfi = 10~^ locates in the region 
^roc > icooi, where Xe is frozen during the isobaric compres- 
sion. Therefore, if the cooling condition tcooi < iff is initially 
satisfied in the post-shock layer, the gas continues to cool 
isobarically until < 100 K, where it produces ~ 10 Mq frag- 
ments. Thus, the boundary of the SMS-forming condition is 
simply given by the requirement tft < tcooi for their initial 
values without the need for considering the recombination 
effect. The dot-dashed line in Fig. |4] presents the condition 
ts = icooi for Xe = 10~^ and in fact coincides with the 
boundary on the low-temperature side (< 8000 K). 



5 CONCLUSION AND DISCUSSION 

In this paper, we have proposed a new scenario for super- 
massive star (SMS) formation in central hot and dense re- 
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gions of the halos formed by the cold accretion shocks in 
the first galaxy formation. Since the gas cools effectively by 
Lya emission in halos with virial temperature Tvir > 10* K, 
location of the accretion shock does not stay at the virial 
radius and shrinks inward. The gas instead flows supersoni- 
cally along cold and dense filaments to the central region 
of the first galaxy, where the flows collide each other to 
produce a hot (> 10'* K) and dense (> 10'' cm~^) ma- 
terial by a shock (Birnboim & Dekel 2003; Keres et al. 
2005; Dekel & Birnboim 2006; Wise & Abel 2007; Greif 
et al. 2008; Wise et al. 2008; Dekel et al. 2009; Bromm & 
Yoshida 2011). We have calculated thermal evolution in such 
a hot and dense region formed by the cold accretion shock. 
For > 8000 K, the efficient Lya cooling allows the post- 
shock gas to cool and to contract isobarically at the value 
of ram pressure from the shock front. To continue the iso- 
baric cooling also below 8000 K, abundant H2 needs to be 
formed and its cooling must be effective. If the density at 
~ 8000 K is high enough (> 10* cm"^) to make the H2 
rovibrational levels to reach the local thermodynamic equi- 
librium, the H2 is dissociated effectively by the collisional 
reaction from excited levels, which suppresses the cooling to 
lower temperature by H2. At this epoch, gravitational in- 
stability of the post-shock layer produces massive fragments 
with > 10^ Mq, which subsequently collapse isothermally 
at ~ 8000 K by the Lya cooling. We have studied ther- 
mal evolution of the post-shock gas for a wide range of ini- 
tial conditions (10^ cm"^ < nn.o < 10^ cm~^ and 3000 
K< To < 10'' K) and have pinned down the conditions lead- 
ing to the isothermal collapse (the hatched region in Fig. [T]): 
To > 6000 (riH.o/lO* cm"^)"* K for nn.o < 10* cm"^ and 
To > 5000 - 6000 K for nH,o > 10* cm-^ for the pre- 
shock ionization degree at a;e,o = 10~^. Since H2 is formed 
by the electron-catalyzed reactions (eq. [HJ, the above condi- 
tion depends somewhat on the initial ionization degree (see 
Fig. [21): for smaller x^fi, the domain of initial conditions 
leading to the isothermal collapse extends towards lower 
density. Those massive clouds continue isothermal collapse 
until very high density ~ 10^® cm""', where they become 
optically thick to the H~ bound-free absorption (Omukai 
2001). The clouds are supposed to collapse directly to SMSs 
without further fragmentation (e.g., Bromm & Loeb 2003; 
Regan & Haehnelt 2009a, b; Shang, Bryan, & Haiman 2010). 
Eventually, the SMSs collapse by the post-Newtonian insta- 
bility, swallowing most of their material, to become seeds of 
SMBHs (Shibata & Shapiro 2002). 

The first galaxies may be enriched with metals to some 
extent as well as dusts dispersed by supernova (SN) explo- 
sions of previous generations of stars. With high enough 
metallicity, the gas can cool to low temperature (~ 100 K) 
by metal- line cooling alone even without H2. In this case, 
the SMS formation would be strongly suppressed. We have 
then repeated the same analysis by considering as well as 
the metal cooling by Cn and Oi. We have found that as 
long as the metallicity is lower than Zcr — 10~^ Zq, the 
metal-line cooling does not change the condition for SMS 
formation from that in the primordial case (see Fig. (3]). Ac- 
cording to some cosmological simulations of the assembly of 
first galaxies (Greif et al. 2010; Wise et al. 2012), the dense 
gas at the center of galaxies is uniformly enriched to ~ 10"'' 
Zq by a pair instability SN of massive population III stars 
with 140 Mq < M < 260 Mq. On the other hand, typical 



population III stars are recently considered to be less mas- 
sive ~ 40 Mq (Hosokawa et al. 2011; Stacy, Greif & Bromm 
2011) and end their lives as ordinary core collapse SNe. In 
this case, the resultant metallicity reduces by a factor of ~10 
(Heger & Woosley 2002; Nomoto et al. 2006) and thus be- 
comes lower than the critical metallicity for SMS formation 
we estimated. 

We here stress that our scenario naturally explains the 
formation of seed BHs in the conditions of first-galaxy for- 
mation without invoking extremely strong UV radiation as 
envisaged in the previous scenario. The necessary condition 
for SMS formation in halos with Tvir > 10* K is the isother- 
mal collapse by the atomic cooling as a result of suppres- 
sion of the H2 cooling in the entire density range. So far, as 
a mechanism to suppress II2 cooling, photodissociation by 
FUV radiation has been considered. In this scenario, how- 
ever, extremely strong FUV intensity J21* > 10^ — 10'' (in 
unit of 10"'^^ erg s~^ cm~^ Hz~^ sr~^) is required to quench 
the H2 cooling (Omukai 2001; Bromm & Loeb 2003; Shang, 
Bryan & Haiman 2010), and halos irradiated by such intense 
radiation are extremely rare (< 10~® at 2; ~ 10; Dijkstra et 
al. 2008). Moreover, if external ionization by cosmic rays 
or X-rays, which promotes the H2 formation, is present as 
well, the FUV intensity needed for SMBH formation is el- 
evated. There would be little possibility (<C 10~'^) for such 
an intense FUV field to be realized in any haloes and thus 
the SMS formation could strongly suppressed (Inayoshi & 
Omukai 2011). On the other hand, in our scenario, collisional 
dissociation, rather than photodissociation, suppresses the 
H2 cooling. Thus, even without FUV radiation, which has 
been considered to be indispensable previously, the isother- 
mal collapse and thus SMS formation can be realized as long 
as the right condition is met for the cold accretion shock. 
Note, however, this mechanism for SMS formation cannot 
operate in all first galaxies since SMBHs are rare objects. If 
we use the number density of halos with mass ~ 10* Mq at 
z ~ 10, ^ 10 Mpc~^ (comoving) and assume each of them 
had a SMS of ^ 10^ the predicted mass density of 

SMBHs ~ 10*^ Mq Mpc"^ (comoving) will exceed the total 
present-day BH mass density ~ 3 x 10^ Mq Mpc~^ esti- 
mated by Yu & Tremaine (2002). Therefore, our conditions 
for SMS formation would be satisfied only in a small frac- 
tion of first galaxies or some other processes, e.g., turbulent 
fragmentation, lack of accreting material etc., suppress this 
mechanism to work for avoiding overproduction of BHs. 

We here briefiy discuss the effect of dust cooling, which 
has not been considered in this paper. If the depletion factor 
of metals to dust grains is as high as the present-day Galac- 
tic value /dep — 0.5, the thermal evolution deviates from 
the isothermal one at > 10^" cm"'' due to the dust cooling, 
if metallicity is higher than Zcr, dust — 10~^ Zq (Omukai, 
Schneider & Haiman 2008). Although this critical metallic- 
ity Zcr, dust is smaller than the critical value due to metal- 
line cooling Zcv — 10"'' Zq by two orders of magnitude, 
the depletion factor in the first-galaxy forming environment 
is highly uncertain. According to theoretical models of dust 
formation and destruction in the first SNe, typically only a 
few % of dust formed at the explosion survives, after be- 
ing swept by the reverse shock, depending on the ambient 
density (Nozawa, Kozasa & Habe 2006; Bianchi & Schneider 
2007). For example, with /dcp — 0.05, the critical metallicity 
becomes Zcr, dust ~ 10 * Zq , which makes the constraint on 



© 0000 RAS, MNRAS 000, 000-000 



8 K. Inayoshi and K. Omukai 



metal pollution less severe. In any case, to predict whether 
the isothermal collapse continues in spite of metal enrich- 
ment, we need more accurate knowledge of the depletion 
factor /dep. 

In this paper, we consider the hot and dense central re- 
gions owing to shocks by the cold accretion flows in the form- 
ing first galaxies. Likewise, a galaxy merging event drives 
inflows, creating a similar environment around the galaxy 
center (e.g., Mayer et al. 2010). If the shocked region satisfies 
our post-shock criterion for the H2 coUisional dissociation, 
the SMS formation is expected also in this case. In merg- 
ing galaxies, however, star formation and also metal enrich- 
ment are expected to have already proceeded significantly. 
Therefore, in the case of the infiows by galaxy merger, SMS 
formation is probably prohibited by the metal-cooling ef- 
fect. The cold accretion shocks in the first galaxy formation 
would provide more easily the suitable conditions for SMS 
formation. Recently, assuming seed BHs with ~ 10^ M!q, 
Di Matteo et al. (2011) and Khandai et al. (2011) discussed 
their growth by the cold accretion flows in the process of the 
flrst galaxy formation. Their results demonstrate that the 
cold flows are less susceptible to feedbacks from the growing 
BHs and high accretion rate is maintained until mass of the 
galaxy reaches > 10^^ where the cold mode of accre- 

tion turns the usual hot virialization mode. As a result, the 
BH is able to grow to > 10^ Mq by 2: ~ 6. Our scenario for 
the SMS formation provides a mechanism for seeding BHs of 
~ 10^ Mq in forming galaxies, which has been assumed in 
studies of Di Matteo et al. (2011) and Khandai et al. (2011), 
while their results complementary demonstrated those seed 
BHs can in fact grow to SMBHs. 

Finally, we remark remaining issues to be explored. In 
our scenario, the physical condition in the post-shock gas 
(especially, its density) is crucial for the SMS formation. 
We have considered a range of density (10^ - lO'^ cm ) and 
temperature (3000 — 10''' K) as the post-shock conditions. 
Currently, we only know the typical values of those param- 
eters (Greif et al. 2008; Wise et al. 2008), but are still lack- 
ing knowledge of the relationship between the post-shock 
conditions and formation conditions (i.e., mass, virialization 
epoch, etc.) of first galaxies. In addition, our assumption 
that the enough mass supply for Jeans instability is avail- 
able through the streaming flow is need further investiga- 
tion. Safranek-Shrader et al. (2010) evaluated the amount 
of accreted gas to be ~ 10^ Mq, which is inhomogeneously 
organized by turbulence (e.g.. Wise & Abel 2007; Greif et 
al. 2008). Therefore, the outcome of the shocked material in 
most of halos could be numerous small sub-regions, rather 
than a massive layer envisaged in this paper. On the other 
hand, even if the turbulent motion dominates, a gravita- 
tional unstable cloud with ~ 10^ Mq could form in the cen- 
ter and collapse subsequently (Wise et al. 2008). As a future 
project, we need to study in what halos the SMS forming 
conditions are satisfied by way of realistic cosmological sim- 
ulations, including e.g., molecular cooling and the radiative 
and chemical feedbacks, for evaluating more quantitatively 
the feasibility of SMS formation in the first galaxies. 

In this paper, we have supposed that, for the SMS for- 
mation, massive clouds must collapse isothermally without 
fragmentation. In fact, three-dimensional hydrodynamical 
simulations by Bromm & Loeb (2003) confirmed in some 
cases that the cloud collapsing isothermally by the atomic 



cooling does not experience fragmentation at least in the 
range < 10® cm~^ and, consequently, a supermassive clump 
forms at the center. Even with some angular momentum, 
fragmentation resulted at most in a binary system in their 
calculation. However, depending on such initial conditions as 
degrees of rotation or turbulence, the clouds would fragment 
into less massive clumps during the isothermal collapse. As 
a future study, it is awaited to clarify the conditions un- 
der which the clouds elude fragmentation by way of three- 
dimensional hydrodynamical calculation. 
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